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ABSTRACT 

We present the results of models of the chemistry, including deuterium, in the 
inner regions of protostellar disks. We find good agreement with recent gas phase 
observations of several (non-deuterated) species. We also compare our results 
with observations of comets and find that in the absence of other processing e.g. 
in the accretion shock at the surface of the disk, or by mixing in the disk, the 
calculated D/H ratios in ices are higher than measured and reflect the D/H ratio 
set in the molecular cloud phase. Our models give quite different abundances 
and molecular distributions to other inner disk models because of the differences 
in physical conditions in the model disk. This emphasizes how changes in the 
assumptions about the density and temperature distribution can radically affect 
the results of chemical models. 

Subject headings: astrochemistry - circumstellar matter — ISM: abundances — 
ISM: molecules — solar system: formation — stars: formation — stars: pre- 
main-sequence 
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Introduction 



Observations and models of the deuterium chemistry in protostellar disks can make 
important contributions to our understanding of the chemistry of the early solar system. 
They can be directly related to the chemistry of the formation of cometary ices, which 
provides a link between the current and primitive solar nebula. Observations of deuterated 
molecules can be used to trace the physical and thermal history of a protostellar disk and, 
combined with models, can determine the relative contribution of the different phases of star 
formation (e.g. molecular cloud, collapse, protostellar disk) to the molecular deuteration 
ratios observed in protostellar disks and in cor nets. Although the elemental abundance of 
deuterium in the Galaxy is only a few xlO~^ (ILinsky et al.l l2006l ) the relative abundance 
of deuterated to non-deuterated molecules can be much higher and depends sensitively on 
the temperature at which the molecules formed, and, at high densities, on the degree of 
molecular depletion. Models of the deuterium chemistry have demonstrated the importance 
of gas phase reactions of deuterated isotopologues of B.^, CHg" and C2H^ in transferring 
deuterium atoms to molecules, as well as the critical role played by grain surface chemistry, 
where reactions of deuterium atoms are able to efficiently form molecules such as deuterated 
formaldehyde and methanol. 



A few deuterated molecules have now been observed in protostellar disks. ICeccarelli et al 



(j2004j ) observed H2D+ in TW Hya and DM Tau and argued that the emission originates in 
the midplane and that therefore H2D"'' provides a means of tracing the ionization level in this 
high densit y, low temperature reg ion where few other molecules exist in the gas phase. A 
later paper (jCeccarelli et al.ll2005l ) reported a detection of HDO in DM Tau where the emis- 
sion is from a region with a temperature of ~ 25 K, far below that at which this molecule 
can be thermally desorbed. This indicates that there is an effic i ent non-thermal desorp - 



tion mechanism acting e.g. photodesorption (IDominik et al.ll2005l : IWillacy &: Langei 



(Note that the observations of both H2D"^ and HD O have been di s puted by 



(|2006|).) DCO+ has been observed in T W Hva flOi et al. 



and in DM Tau (IGuilloteau et al.ll2006l ). Ivan Dishoeck et al. 



2008 



0.035. 



Oi et al 



2000). 



Guilloteau et al. 

2003h 



( 20031 



van Dishoeck et al 

T 

used JCMT and found 



a beam-averaged ratio of DCO'^/HCO"*^ 
a more complex distribution of DCO+ which does not follow that of HCO" 



((20081) using the SMA found 
In this case 

DCO+/HCO+ increases with radius between 30 and 70 AU from 0.01 to 0.1, and DCO+ 
disappea rs rapidly at -R > 90 AU . In contrast HCO"^ and CO are present out to 20 AU. In 
DM T au [Guilloteau et all J2006h find DCO+/HCO+ ~ 4 xlO^^ in the outer disk, loi et al. 



(120081 ) also observe DON in TW Hya and determine DCN/HCN to be 0.017. 



These observations are all for the outer regions of the disk (i? > 30 - 50 AU). Here our 
focus is on i? < 30 AU. The inner regions of disks are of great interest because they cover 



3 



the region where planet formation took place in the solar nebula. They are now becoming 
acces s ible to observational investigation e.g. with the Spitzer Space Telesco pe (ICarr fc Naiita 



2008 : ISalvk et all l2008l : iLahuis et all 120061 ) or the Keck Interferometer flGibb et al 
20041 ). In the future, ALMA will further open up this region for study. 



2007 



solar system data provides an additional constraint on the models. Observations of 
deuerated molecules in comets are limited to two molecules (HDO and DON) in four comets 
- Halley, Hale-Bopp, Hyakutake and C/2002 T7 LINEAR. The ratios are sirnilar in all of the 
observed comets with HDO/H2O ~ 5 - 7 x 10~ '^ ( Bockelee-Morvan et al .11 1998 ; Eberhardt et al 



1995 



Balsiger et aDll995l : Elutsemekers et al.ll2008h and DCN/HCN ~ 2 xlO'^ dMeier et al. 



19981 ). These values are higher than the solar elemental D/H ratio. The similarities between 
molecular D/H ratios observed in the interstellar medium and those seen in comets have led 
to the suggestion that cometary ices have their origin in the interstellar medium. However, 
the re is sufficient differ ence between the two data sets for this to still be an open question 



e.g 



Bergin et al.l (120071 ) and to suggest that processing during star formation, or in the early 
solar nebula could affect the ratios. Once comets formed they are thought to have undergone 
relatively little subsequent processing and therefore their abundances reflect the composition 
of material at the time of their formation. 

Enhanced D/H ratios are seen elsewhere in the solar system. For examp le, interplanetary 



dust particles are enri ched in D with respect to the terrestrial value (e.g. [Messenger et al. 



2003 



Messenger! |2000| ) . They show an extremely wide range of D/H ratios fiRobert et al. 
2OOOI ) but the carriers have not yet been identified. Enhancements are also seen in primitive 
meteorites (LL3 and carbonaceous chondrites, as well as meteoritic water), but the anomalies 
are generally smaller. The observations do not yet allow for th e origin of the deuterated 
material to be defini tively identified, and both interstellar (e.g. iPizzarello fc HuangI |2005| ) 
and protostellar (e.g. iRemusat et al.ll2006l ) origins have been claimed. 



The deuterium content of the giant planets Jupiter and Saturn is considered to be a 
relic of the early protosolar nebula, since deuterium is neither formed nor destroyed in these 



planets, and 


their 


(Owen et al. 


1986) 



mation about the deuteration levels in the protosolar nebula at their time of formation. 
Measurements in Jupiter suggest t hat the D/H ra t io is e nhanced over the interstellar value 



elemental ratio by a factor of a few. iMahaffy et al.l (119981 ) used the Galileo Probe Mass Spec- 
trometer to d erive (D/H.)^^, = 2.6 ± 0.7 x 10~^, in close agreement with the ISO-LWS value 



of 2.2 X 10~^ ( iLellouch et al.ll200ll ). Both these values depend on the models used to analyze 
the data. A direct measu r emen t of D/H in the limb of Jupiter using Lyman-a emission was 
made by lBen Jaffel et al.l (119981 ) who found (D/H)//2 = 5.9 x 10~^. Deuterated methane has 
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also been observed wi th a ratio (D/}i)cH. = 1-8 xlO ^-2.9x10 ^ (iFeuchtgru 
Encrenaz et ali 1999 ). Saturn has similar D/H ratios to Jupiter (1.7 x 



ber et ah 1999 



0" 



Lellouch et a 



Bezard et al 



19971 ) 



(|2001| )). In N eptune and Uranu s the ratios are higher - 6.5 x 10 ^ 
or 5.6 xlO~^ (lOrton et aLlll992l ) for Neptune and 5.5 x 10~^ (IFeuchtgruber et al.lll999l ) in 



Uranus, possibly as a consequence of the inclusion of highly deuterated icy planetesimals 
into the planets as they formed. 



In our previous paper (jWillacyi 120071 ) we considered the deuterium chemistry in the 
outer regions {R > 50 AU) of a protostellar disk. We showed that although the gas phase 
molecular D/H ratios can be considerably altered by chemical processing in the disk, the 
grain mantle ratios reflect those set in the parent molecular cloud, supporting the idea that 
comets may contain interstellar material that survived the formation of the Sun and planets 
relatively unchanged. Here we extend our previous work to cover the inner disk from 0.5 - 30 
AU. This radius range covers the comet and planet forming region, and is also appropriate 
for comparison to infrared observations of the inner regions of protostellar disks. 

We present the results for a typical T Tauri star disk and for a more massive disk 
with a higher surface density similar to the minimum mass solar nebula (MMSN). Our 
models incl ude multiply deuterated molecules since in high density regions these can become 
important fICeccarelli fc DominikI l2005l : [Roberts et al.l 120021 . |2003| ) although their effects in 
the warm inner disk are expected to be less than in the colder outer disk. We compare our 
results with the available observations of the inner disk (non-deuterated molecules) and with 
observations of deuterated molecules in comets. 



The model 



To fully model a disk self-consistently requires many chemical and physical processes 
to be taken into account. To simplify this process and to reduce the required computing 
time we separate the calculation of the physical structure of the disk from the chemical 
network. We obtain the disk structure (dens i ty an d grain temperature) from a pre-existing 
hydrodynamical model (Id'Alessio et al.lll999l . l200ll ). These models use the a prescription to 
determine the density and temperature distribution throughout the disk. 

We consider two of d'Alessio's models with different mass accretion rates. The first 
(Model 1) has M = 10"^ Mq yr'^ a central star of mass = 0.7 Mq, radius R* = 2.5 R© 
luminosity L,, = 0.9 Lq and temperature T* = 4000K. The viscosity parameter a is 0.01. 
The dust distribution is described by a power law n{a) = noa~^ where a is the grain radius 
and p = 3.5. The maximum grain size is 0.25 /im and the dust is assumed to be well mixed 
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wi th the g a s. Th e surface density at 5 AU is 38.8 g cm ^. This model is the same as used 



m 



Willacyl (120071 ) (hereafter W07) and so represents a continuation of the model to smaller 



radii. 

The mass of the disk in this model is somewhat lower than that expected in the early 
solar system. We therefore also consider a more massive dis k (Model 2) w here the surface 
density is closer to that of the minimum mass solar nebula (jHayashilll98ll ). In this model 
M = 10~® Mq yr~^ and a = 0.025. All other parameters remain the same as in Model 1. 
The surface density at 5 AU for Model 2 is 150 gcm~^. 

Although the hydrodynamical models assume that the ga s and grain temperature ar e 
equal, this is not th e case in the surface layers of the disk e.g. iKamp fc DuUemondl (120041 ): 
Gorti fc HoUenbachl (120041 ). where the gas temperature can be much higher than the grain 
temperature. Since the gas temperature can critically affect the chemical (and especially 
the isotopic) abundances we el ect to calculate the gas tern perature separately. This calcu- 
lat ion is based on th e work of iKamp fc DuUemondl (120041 ) and is described in more detail 



m iWoods fc Willacyl (120091 ). We summarize here the processes that are included. Heat- 
ing is provided by the photoelectric effect, photodissociation of H2, the formation of H2 on 
the surfaces of dust grains, the ionization of carbon atoms, cosmic rays and stellar X-rays. 
Gas-grain collisions can also act as heating processes in regions where the grains are hotter 
than the gas, but this is not an important process in the inner disk, where this process is 
more likely to cool the gas. Cooling occurs via the emission lines of atomic oxygen and 
ca rbon, C"*", CO and CH, a nd by Lyman-a cooling. The rates for these processes are given 
in 



Woods fc Willacvl (I2009D . 



For H2 formation heating we use the H2 formation rate from lCazaux fc TielensI (120021 ). 
which assumes that H atoms can both chemisorb and physisorb to bare grains. Chemisorp- 
tion only occurs on bare grains and this rate is therefore only valid at the top of our disk. 
Below this where the grains are covered in ices, we assume that only physisorption occurs 



with the binding energy given in ICazaux fc Tie- 



end (I2OO21) . This results in cooler tempera- 



tures than were found by iKamp fc DuUemondl (120041 ) in regions just below the disk surface 
where ices have begun to form. Figure [T] shows the density and temperature distributions 
for the two models used. 

The final parameter re quired is the UV field. In order to estimate the UV field in the 



disk we adopt the model of iRichling fc Yorkd (120001 ) which uses a ray tracing method and 
the fiux limited diffusion approximation to determine the field at each point in the disk This 
allows us to model the effects of UV irradiation by the interstellar UV field, and by the central 
protostar and to include scattering of t he UV photon s by d ust g rains. The graiii absor ption 
and scattering coefficients are based on lDraine Sz Led (119841 ) and iPreibisch et al.l (119931 ) who 
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Fig. 1. — The physical parameters used in (a) Model 1 (a = 0.01) and (b) Model 2 (a = 
0.025). Left: The gas density; Center: the grain temperature; Right: the gas temperature. 
Gas density and grain temperature are supplied by the models of d'Alessio, and the gas 
temperature is calculated separately by balancing the heating and cooling processes. 



Density Grain temperature Gas temperoture 
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assumed a mixture of silicates, amorphous carbon and dirty ice-coated silicate grains. The 
disk is divided into square cells (Az = AR = 0.02 AU) and the radiative transfer equation 
is solved separately for both s ources of radiation and for the direct and diffuse (scattered) 
components of each field (see iRichling fc Yorkd l2000l . for details). The total UV flux at a 
given point is the sum of the contributions from the stellar and interstellar fields. In general, 
this results in a total UV field that is higher than that esti mated from our previous slab 
method: see for example W07. A similar effect was found by Ivan Zadelhoff et al.l (120031 ) for 
a more complex radiative transfer model in the outer disk. 

Using the density and temperature distributions provided by the d'Alessio disk model 
we calculate the UV field for both the interstellar and stellar fields. For the latter we assume 
that the UV field emitted by the T Tauri star is 500 Go at 100 AU where Gq is the standard 
interstel lar radiation field (ISRF). This field strength is based on observations of T Tauri 
stars by lBergin et al.l (120031 ). but we do not take into account the possibility that the T Tauri 
spectrum has a different shape to the ISRF and that it may be dominated by strong emission 
features such as Lya. This can significantly impact the chemistry of the surface layers of 
the disk sinc e Lya can dissociate some molecules such as OH and CH4, but not others. 
Bergin et al. J2003h found that this effect can accou nt for high CN/HCN ratios observed in 
some disks (IDutrey et al.lll997l : iKastner et al.l 119971 ). 



Following previous authors (e.g. lWoods fc Willacyll2009l : lMarkwick et al.ll2002l : IWillacy et al 



19981 ) we assume that material advects inwards towards the star. We do not have a com- 
pletely integrated dynamical and chemical model and instead follow parcels of gas, starting 
at the outer edge of the disk (in this case 35 AU) which gradually move inwards. The parcels 
move inwards in discrete radial steps of 0.5 AU and travel along lines of constant gas pressure 
scale height. The time that a gas parcel remains at a given grid position can be determined 
from the radial velocity: 

M 



V,. 



2nT.R 

where S is the surface density of the disk at radius, R. The timescale for accretion. At 
AR/vr and therefore the time that a parcel of gas spends at each radius is given by 



'1 



At 



2txT.RAR 
M 



(2) 



S scales roughly with 1/i? for i? > 10 AU (Id'Alessio et al.l I2OOII ) so equation [2] has only 
a small dependence on R. Starting at 35 AU, a parcel of gas will migrate into the star in 
approximately 0.41 million years. In our grid, the parcel spends roughly 6000 years at each 
radius before moving inwards to the next gridpoint. We have 70 radial and 56 vertical grid 
positions. 
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2.1. The chemical network 



Our chemical network is the same as used by WO 7. We t ake the basic gas phase reaction 
set from the UMIST database RATE95 ( Millar et al. 1998h and add deuterated reactions 
using the techniques of iMillar et al.l (119891 ) and lRodgers fc Millarl (jl996l ). The reaction rates 
for th e multiply deuterated isotopologues of Hg 



and CHj]" are taken from [Roberts et al. 



( l2004j ). For details of the grain surface chemistry see W07. The reaction network links 
227 gas phase species (of which 115 are deuterated) and 91 grain species (44 of which are 
deuterated) by 9489 reactions. 

We also include freezeout, thermal desorption and cosmic ray heating of grains. The 
rate at which a species, X, will freezeout onto dust grains is given by 

(3) 



Hg is the 



k freeze = SxCx < Tia Ug > Vx S ^ 

where Sx is the sticking coefficient (= 0.3 for all species), a is the grain radius 
number density of grains and vx is the gas phase velocity of X. Cx is a factor to take into 
accou nt the increase in freezeout r ate for positive ions encountering a negatively charged 
grain (lUmebayashi fc Nakandll980l ). Cx is 1 for neutral species, and 1 + 16.71 x 10~^/ (aT^^) 
for positive ions, where Tg^ is the grain temperature. Ions are assumed to recombine on the 
grain surface in the same way as they do when reacting with electrons in the gas. All species 
are assumed to freezeout with the exception of He and He+ which have a very low binding 
energy and are easily thermally desorbed even at low temperatures. Any He"'' that hits a 
grain is immediately returned to the gas as neutral He. 



The thermal desorption rate is given by 



k 



therm 



Z/qC 



(4) 



z/q is the frequency of oscillation between the absorbate and the grain surface and is given 
by 

= ^/l2nsED/n'^m) (5) 

where Ug is the surface density of sites (~ 1.5 xlO^^ cm~^), m is the mass of the accreting 
species and E^, is its binding energy. The binding energies used are summarized in Table [T] 
For deuterated species, the binding energies are assumed to be 21 K more than their un- 
deuterated equivalents (21 K is the differ ence in zero point energy between the hydrogen 
and deuterium atoms) JCaseUi et alll2002l ). For example, ^^(HaO) = 5770 K, Ed(HDO) = 
5791 K and -5'n(D?0) = 5812 K . Bin ding energies for specie s not listed in Tabled] are taken 
from iHasegawa fc HerbstI (119931 ) and iHasegawa et al.1 (119921 ). 



Cosmic ray heating can remove weakly bound molecules from grains in regions where 
the temperatures are low, and where the surface density of the disk is low enough that 
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cosm ic rays can penetrate. Th e desorption rates due to cosmic ray heating are taken 
from iHasegawa fc HerbstI (Il993l ) using updated binding energies for some species e.g. CO 
dOberg et al.lboosh . 

In our outer disk models (W07) we included photodesorption since at i? > 50 AU the 
grain temperatures are low enough that thermal desorption is not efficient for most species. 
Photodesorption was found to remove molecules in the surface layers of the disk even when 
the grain temperature was low. In the inner disk model the grain temperatures in the surface 
layers are warm enough that even the molecules with the highest binding energies can be 
removed efficiently by thermal desorption and therefore photodesorption will have little effect 
on the molecular abundances. Hence we have chosen to exclude it from our current models. 



2.1.1. Ionization processes 



Ionization in the disk can arise from several processes including irradiation by UV 
photons (both interstellar and stellar), protostellar X-rays and cosmic rays. There is also a 
contribution from the decay of radioactive isotopes. 

The UV field at each disk posit ion is calculated a s described in Section [21 The stellar 
field is taken to be 500 Gq at 100 AU (IBergin et al.ll2003l ) where Gq is the standa rd interstellar 



radiation field. We assume that the stellar field does not dissociate CO and H2 jAikawa et al. 



2OO2I ). although these molecules are dissociated by the interstell a r field . We calculate the 
self-shielding of CO and H2 using the slab model of iLee et al.l (119961 ). This model was 
developed for molecular clouds and provides data for the shielding due to H2, CO and dust as 
a function of column density. The linewidths assumed are 3 kms~^ which is much larger than 
observed in disks, where the velocity dispersion in the outer disk at least, is almost thermal 
(IGuillqteau fc Dutreyl Il998l ). To take account of this we have followed lAikawa fc Herbst 



f ll999bl ) in scaling the column densities in Table 10 of iLee et al.l by Cj,/3 kms ^ where Cg is 



the local sound speed in the disk. This scaling factor is only required for H2 since the CO 
dissociation lines are broader due to predissociation. 

HD does not self-shield, but some of its lines do o verlap with those of H2. This means 
that the H2 can provide some shielding. iBarsuhnI (119771 ) estimated this would reduce the HD 
photodissociation rate by 1/3, assuming that the overlapped HD lines are totally shielded 
by H2. Here we modify the HD dissociation rate accordingly. 

Cosmic rays can also cause ionization, both directly and indirectly (by interaction with 
H2 producing UV photons). The rates for both processes are taken from the UMIST ratefile, 
with an assumed cosmic ray ionization rate Co = 1-3 x 10~^^ s~^. Cosmic rays can produce 
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Table 1. The binding energies (Ed) used to determine the thermal desorption rates of the 

abundant mantle species. 



Species 


Ed (K) 


Reference 


Species 


Ed (K) 


Reference 


H 


600 


1 


D 


621 


2 


H2 


315 


8 


C 


800 


3 


CH 


645 


9 


CH2 


956 


9 


CH3 


1158 


9 


CO 


855 


4 


CO2 


2860 


5 


H2C0 


1760 


3 


CH3OH 


4240 


7 





800 


3 


O2 


1210 


3 


OH 


1259 


9 


H2O 


5770 


6 


N 


800 


3 


N2 


790 


4 


NH 


604 


9 


NH2 


856 


9 


NH3 


3080 


5 
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ionization in the disk only if the surface density (S) is lo w enough for them to penetrate i.e. 
S is less than 150 g cm~^ (jUmebayashi fc Nakandll98ll ). We assume that cosmic rays only 
penetrate the disk vertically, and that they can penetrate from both above and below the 
disk. Their ionization rate is given by 



CcR = -Co[exp(-Si(z,i?)/102) +exp(-S2(z,i?)/102)] 



(6) 



(jSemenov et al.ll2004l ) where Si is the surface density between the height above the midplane, 
z, and the top of the disk, and S2 is the surface density between z and the bottom of the 
disk. 

The decay of radioactive nucleides is an additional source of ionization. ^^Al can de- 
cay to form excited ^^Mg, which in turn decays by either positron decay or by electron 
capture. The ionization rate of these processes is assumed to be Cai = 6.1 x 
dUmebavashi fc NakanolllQSlh . 



10-19 s-1 



A final source of ionization are X-rays of which T Tauri stars are strong emitters. We 
use the approach of iGorti fc HoUenbachl (120041 ) to calculate the X-ray photoionization rates 
of atoms and molecules. This method assumes that the X-ray photoionization rate of a 
molecule is the sum of the rate for its constituent atoms and that X-ray ionization leads to 
the loss of a single electron. The ionization rate per atom depends on the incident X-ray 
flux and the cross-sect ion of the atom. We calc ulate fits to the energy dependent atomic 
cross-sections given in I Verner fc YakovlevI (119951 ). The attenuation depth of X-rays is very 
small and X-ray ionization is most efficient in the surface layers of the disk. For further 
details, see IWoods fc Willacvl J2009h . 



2.2. Recombination of N2H"'" 

The recombination of N2H+ with electrons has recently been the subject of study in the 
laboratory with two different groups producing two different results. There are two possible 
pathways for this reaction 

N2H+ + e- — ^ NH + N (7) 
— ^ N2 + H (8) 

The pathway followed by this reaction can influence the abundance of nitrogen-bearing 
molecules in the gas phase. If the primary product is NH then in cold regions of the disk, 
NH will quickly freezeout and be converted to NH3 ice, removing nitrogen from the gas 
phase. N2, on the other hand, does not react on the grains and is relatively easy to desorb. 
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N2 is also an important means of controlling the degree of ionization fjCeccarelli fc Dominik 
20051 ). The presence of N2 in the gas can prevent the transfer of deuteration along the chain 
of isotopologues from H3 to D3 by destroying the less deuterated isotopologues before they 
have a chance to form Dt. 



In W07 we took the results of iGeppert et al.l (120041 ) which found that the NH route is 
most likely to happen, o ccurring 65% of the time. However, more recent laboratory work 
from iMolek et al.l (120071 ) has cast some doubt on this result. iMolek et al.l (120071 ) find that 
the route leading to N2 dominates, with no significant NH being formed. 



Using the iGeppert et al.l result s leads to lower abundances of N2 in the midplane com- 
pared to models using iMolek et al.l This will translate into a higher deuteration of Hjj' and 
lower abundances of gas phase nitrogen molecules. 



Here we use the more recent results of iMolek et al.l and assume that little NH is formed 
by the recombination of N2H"'". 



2.3. The initial conditions 



We assume that the material incorporated into the disk has first been processed to some 
extent in the parent molecular cloud. The input abundances for the disk model are therefore 
taken from the output of a molecular cloud model. The latter uses the elemental abundances 
found in Table [2] and is allowed to run for 1 Myrs at a total hydrogen density uh = 2 xlO^ 
cm~^, temperature = 10 K and visual extinction Ay = 10 magnitudes. All the chemical 
processes included in the disk model are also included in the molecular cloud model with the 
exception of ionization by the decay of radioactive nuclei. The input disk abundances are 
given in Table [31 The processes that drive deuteration ch e mistry in molecular clouds hav e 
been well described by previous authors e.g. iRoberts et al.l (120041) ; iRoberts fc Millarl (I2OOOI ) ; 



Charnlev et all (Il997f ): 



Brown fc Millai] (Il989l ) 



Tielend (119831 ) 
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Table 2. Elemental abundances used in the model. Abundances are given with respect to 
the total abundance of hydrogen atoms, nn = 2 n(H2) + n(H). Initially we assume that 1% 
of the hydrogen is atomic, with the rest molecular, and that 1% of the deuterium is atomic, 

with the rest in HD. 



Element 


Abundance 




n{x)/nYi 


H2 


0.495 


H 


0.01 


He 


0.14 


D 


1.6 X 10-^ 





1.76 X 10"^ 


C 


7.3 X 10^5 


N 


2.14 xlO-5 


Fe 


3.0 xlO-9 


Mg 


7.0 xlO-9 
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Table 3. Input fractional abundances with respect to total number of hydrogen nuclei for 
the disk model as determined by a molecular cloud model at 1 Myrs. The results for two 
models are shown. The numbe rs vary a little from those in Table 3 in WO 7 because of th e 



use of the iMolek et al.l (120071 ) results for N2H"'" recombination instead of iGeppert et al. 

( l2004j ). The physical parameters used were density, uh = 2 x 10^ cm~^, temperature = 10 
K, and Ay = 10 magnitudes. Freezeout, thermal desorption and cosmic ray heating of 
grains were included. Grain chemistry contributes to the formation of some of the ice 

molecules e.g. water, whereas others e.g. CO form in the gas phase and then freezeout but 

are not produced in situ on the grain surfaces. 



Molecule 



Abundance 



Molecule 



Abundance 





Gas 


Grain 




Gas 


Grain 




2.8 


(-5) 






6.7 


(-7) 




H3+ 


3.0 


(-9) 




H2D+ 


5.8 


(-10) 




HD+ 


9.8 


(-11) 






1.5 


(-11) 




CO 


3.2 


(-5) 


2.7 (-6) 


CO2 


5.8 


(-8) 


2.2 (-7) 


HCO+ 


3.5 


(-9) 




DCO+ 


3.1 


(-10) 




H2CO 


2.4 


(-7) 


1.6 (-6) 


HDCO 


1.6 


(-8) 


7.7 (-8) 


CH3OH 


1.6 


(-10) 


4.0 (-8) 


CH3OD 


9.3 


(-12) 


9.1 (-10) 


CH2DOH 


1.5 


(-11) 


3.2 (-9) 





8.5 


(-7) 




O2 


4.4 


(-8) 




OH 


2.4 


(-8) 




OD 


4.1 


(-8) 




H2O 


3.8 


(-8) 


1.4 (-4) 


HDO 


3.8 


(-9) 


1.2 (-6) 


D2O 


4.2 


(-11) 


1.3 (-8) 


N 


1.0 


(-7) 




N2 


1.9 


(-6) 


6.5 (-8) 


NO 


4.5 


(-8) 




CN 


6.6 


(-8) 




HCN 


2.1 


(-8) 


2.4 (-6) 


DON 


6.7 


(-10) 


2.0 (-8) 


HNC 


9.9 


(-9) 


1.1 (-7) 


DNC 


3.9 


(-10) 


2.1 (-9) 


HC3N 


2.5 


(-7) 


9.9 (-7) 


DC3N 


9.4 


(-9) 


3.3 (-8) 


NH3 


8.6 


(-9) 


1.3 (-5) 


NH2D 


2.7 


(-10) 


6.7 (-8) 


NHD2 


5.1 


(-12) 


3.1 (-10) 


ND3 


8.7 


(-13) 


2.1 (-11) 


N2H+ 


2.8 


(-10) 




N2D+ 


2.8 


(-11) 




C 


3.9 


(-7) 




C+ 


3.1 


(-9) 




C2H 


6.4 


(-9) 




C2H2 


1.1 


(-7) 


3.3 (-7) 


CH4 


3.1 


(-6) 


1.1 (-5) 


CH3D 


1.6 


(-7) 


3.4 (-7) 


CH2D2 


6.1 


(-9) 


1.8 (-8) 


CHD3 


1.5 


(-10) 


5.0 (-10) 


CD4 


8.2 


(-13) 


2.1 (-12) 


HD 


1.2 


(-5) 
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Results 



In common with previ o us authors e.g. lAikawa fc HerbstI (jl999bl ): IWillacy fc Langer 



((20001); lAikawa et all (120021 ): Ivan Zadelhoff et all (120031 ) we find that the disk can be di- 



vided vertically into three chemically distinct layers. Grains are coldest in the midplane, 
leading to the existence of ices of at least some species at i? > 2 AU. Inside of 2 AU the 
grains are hot enough that all molecules are desorbed in the midplane. Above the midplane 
is a molecular layer where thermal desorption is efficient and the disk is sufficiently optically 
thick that molecules can survive in the gas phase. This leads to high gas phase molecular 
abundances. In the surface layer UV and X-rays can penetrate, dissociating molecules into 
their constituent atoms and ions. The distinction between the midplane and the molecular 
layer is not as clear as that seen in outer disk models, since the grains are warm enough even 
in the midplane for some volatile species (e.g. CO and N2) to desorb. We use these three 
layers in the discussions below. 

Figure [2] shows the fractional abundance distributions in the inner disk. The abundance 
distribution is a combination of the effects of the formation and destruction processes at a 
given radius, together with the transport to that radius of molecules formed further out in 
the disk. Many molecules desorb as their sublimation temperature is reached and then show 
little variation in their gas phase abundance with radius e.g. N2 and CO show no significant 
change in abundance across the inner 30 AU of the disk below the surface photodissociation 
layer. The upper extent of these molecules is determined by photodissociation. N2 shows a 
more complicated behavior than CO, with a low abundance region (where x(N2) < 10~^^). 
Here photodissociation is still efficient and the resulting nitrogen atoms are incorporated 
into other nitrogen-bearing molecules e.g. HCN. HCN has a high binding energy of 4173 K 
and the grain temperature in this region is low enough that once HCN freezes out it is not 
able to thermally desorb. This leads to a loss of nitrogen atoms from the gas and hence 
N2 cannot reform. Above the low abundance layer, HCN still forms from nitrogen atoms, 
but the difference is that the grains are warm enough for desorption of HCN to be efficient. 
Hence nitrogen atoms can cycle between N2 and HCN. Below the low abundance region N2 
does not photodissociate and has a high abundance. 

The effects of an increase in temperature towards the disk surface and towards the star 
can be seen in the increase in abundance of less volatile species such as H2O and HCN in 
these regions. Ammonia has a lower binding energy and is present for a larger radial extent 
than either H2O or HCN, but it does not survive in the surface layers where it is quickly 
destroyed by reaction with abundant CN, producing HCN. The HCN freezes out and is 
retained on the grains and NH3 is not reformed. NH3 only exists in regions shielded from 
UV radiation. 
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Fig. 2. — Fractional abundances in the inner disk for Model 1 (a = 0.01). Note change of 
radial scale for molecules that are only abundant in the inner few AU i.e. H2O, HDO, NH3 
and NH2D. The dashed line indicates the surface of the disk (at 6 pressure scale heights 
above the midplane). Also marked are optical depths, r, of 0.1, 1 and 10. 
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H2CO shows two abundance peaks in the mi dp lane, as well as a peak in the molecular 
layer above it. The outermost midplane peak (at R = 20-23 AU) and the molecular layer 
peak are due to thermal desorption. Between 7 and 20 AU H2CO is quickly destroyed by 
gas phase reactions with HCO''' and DCO"*" followed by neutralization of the resulting ions 
(HaCO^ and H2DC0^) either with electrons or on grain surfaces. While some of the ions 
will reform H2CO, 1/3 will instead form CO, leading to a removal of H2CO from the gas. 

The second peak at i? < 8 AU is due to formation in the gas phase by 

O + CH3 — > H2CO + H (9) 

The oxygen atoms are supplied by the destruction of CO by reaction with He"*". The CH3 is 
supplied by the breakdown of C3H4 (which is desorbed at ~ 7 AU). HDCO forms by similar 
processes. 

C2H exists in the boundary between the molecular and surface layers as a photodisso- 
ciation product of C4H2 at i? > 10 AU and of HC3N inside of this radius. 

The abundances of molecular ices do not vary much across the disk, until their desorption 
temperature is reached. The ice abundances are set by reactions in the molecular cloud and 
are not altered by grain chemistry in the disk. 

3.1. Gas phase column densities 

Figure [3] shows the radial variation in column density measured vertically from the 
midplane for some molecules as calculated in Model 1. Most molecules show a sharp increase 
in column density due to thermal desorption as the radius decreases e.g. HCN, CO2, NH3. A 
corresponding increase in fractional abundance can be seen for these molecules in Figure [2J 
The increase in the column density of the equivalent deuterated molecules occur at slightly 
smaller radii because of their slightly higher binding energies. The fall in DCN at i? < 
1.5 AU is because of its reaction with hydrogen atoms to form HCN. The reverse reaction 
(D + HCN — s- DCN + H) can also occur, and at the same rate, but DCN does not reform 
because the atomic D/H ratio is low at small radii. iV(CN) also increases closer to the star. 
It is a photodissociation product of DCN and HCN and therefore exists mainly in the surface 
layers. 

H2CO shows two column density peaks: one around 23 AU and one around 3-4 AU 
corresponding to the fractional abundance peaks in the midplane seen in Figure [21 
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Fig. 3. — The calculated column densities in Model 1 of some molecules in the inner 30 AU. 
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Table 4. Calculated vertical column densities from midplane to the disk surface for some 

gas phase species in Models 1 and 2. 



Molecule 




Model 1 ( 


a = 


0.01) 




Model 2 (a = 


0.025) 




1 AU 


5 AU 


10 AU 


1 AU 


5 AU 


10 AU 


H2 


1.6 


(25) 


2.4 


(24) 


2.3 


(24) 


4.9 (25) 


1.9 (25) 


9.2 


^24) 


CO 


1.1 


(21) 


3.8 


(20) 


3.6 


(20) 


1.9 (21) 


1.9 (21) 


6.0 


^20) 


CO2 


3.5 


(19) 


1.2 


(18) 


7.0 


(14) 


7.3 (19) 


1.7 (19) 


6.8 


:i3) 


H2O 


4.1 


(21) 


8.2 


(16) 


2.5 


(16) 


1.4 (22) 


1.4 (17) 


7.1 


:i5) 


HDO 


4.2 


(19) 


6.4 


(12) 


8.7 


(11) 


1.7 (20) 


8.7 (12) 


2.6 


:ii) 


OH 


2.4 


(17) 


7.6 


(16) 


6.8 


(16) 


1.2 (17) 


1.1 (17) 


8.3 


'16) 


HCN 


1.1 


(20) 


6.6 


(14) 


4.2 


(14) 


4.1 (20) 


1.8 (14) 


3.6 


:i4) 


DCN 


5.0 


(14) 


2.4 


(11) 


4.3 


(10) 


2.6 (15) 


1.4 (10) 


2.6 


:io) 


CN 


1.5 


(18) 


2.2 


(17) 


8.6 


(16) 


3.4 (18) 


1.9 (17) 


8.3 


'16) 


C2H2 


2.0 


(19) 


6.0 


(18) 


3.6 


(16) 


3.3 (19) 


2.2 (19) 


2.5 


:i7) 


CH4 


6.5 


(18) 


2.3 


(18) 


4.5 


(18) 


2.8 (18) 


4.8 (18) 


9.5 


:i7) 
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3.2. Deuteration in the inner disk 



In low temperature regions such as molecular clouds enhancements in molecular deuter- 
ation are driven by the reactions of H2D"'", 01320"*" and C2HD. These form by 



CH+ 
C2H+ 



HD ^ H2D+ + H2 + AE1 
HD ^ CH2D+ + H2 + AE2 
HD ^ C2HD+ + H2 + AE3 



(10) 
(11) 



where AEi, AE2 and AE3 are activation barriers, which inhibit the reverse reactions at low 
temperatures. {AEi = 220 K at low temperatures, falling to ~ 130 K as the temperature 
increases, AE2 ~ 370 K and AE^ ~ 530 K. This allows deuteration to be transferred to 



other molecules via ion-molecule reactions. [Roberts et all (120031 ) and [Roberts et al.l (120041 ) 
showed that multiply deuterated forms of H3 , CH3 and C2H^ also play an important role 
in determining the deuteration i n cold r e gions . Grain surface re a ctions can also efficiently 
increase deuteration in ices e.g. iTielenj (119831 ): iBrown fc Millarl (119891 ). For example, the 
high abundances o f deuterated formaldehyde and methanol ob served in star forming regions 



(jParise et al.l 120021 : iLoinard et al.l l2002l : iBacmann et al.l |2003| ) have been attributed to the 



reaction of deuterium atoms on grains during an earlier low temperature phase of evolution. 

In our disk model we find that the isotopologues of and CH;^ are the main drivers 
of deuteration with some contribution from the reaction of deuterium atoms. 02110"*" plays 
only a minor role. Figure H] shows the fractional abundances of and its isotopologues, 
and the level to which this molecular ion is deuterated is shown in Figure [51 The abundance 
of Hj]" peaks above the midplane where cosmic ray ionization is efficient and hydrogen is 
mainly atomic. But here the temperatures are relatively high, and the deuteration is low. 
In the midplane, the D/H ratios fall off rapidly with decreasing radius, and the more highly 
deuterated forms are present only at larger radii (Fig. [5]) . In the outer disk [R > 100 



AU), the models of W07 and ICeccarelli fc DominikI (120051 ) found that HD^ and Dj]" are 



produced with significant abundances and therefore play an important role in the deuterium 
chemistry of this region. In the inner disk the temperature is higher and the abundances of 
these isotopologues are much lower. Deuteration is therefore driven only by H2D"*". 

In common with Hj]", the deuteration of CHj]" and atomic hydrogen peak in the midplane 
at R > 15 AU. The abundance of CHj]" is very small in this region, however, and its iso- 
topologues play little part in driving the deuteration. Deuteration by reaction with atomic 
D is most efficient ai R > 17 AU, and inside of 5 AU. The abundance of atomic D and H 
in the midplane is relatively high because of the high grain temperature. This reduces the 
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residence time of the atoms on the grains so that they are more hkely to desorb than to react 
with other atoms or radicals. 

Other species do not necessarily show a decrease in D/H ratio with decreasing radius 
(and increasing temperature). Many molecules retain the D/H ratio that they had when 
they formed in the cold molecular cloud or in the colder regions of the disk, as they are 
transported inwards. The transport process occurs on a shorter timescale than the chemical 
reactions that would destroy the deuteration. For example, the deuteration of H2CO and 
NH3 reflects that set in the ices during the molecular cloud phase and is not affected by 
chemistry in the disk. 

The deuteration of gas phase water is less than that seen in the ice. Water is desorbed 
in two regions, in the surface layers and at i? < 2 AU. At high z HDO is formed by 

OD + H2 — >HDO + H (12) 

and destroyed by 

HDO + H — ^OH + HD (13) 

Similar reactions affect H2O, but since the destruction product is OH, this is able to cycle 
back to form more water by reaction with H2. Hence a decrease in HDO/H2O is seen at high 
z. Below this the ratio is higher (yellow stripe in Fig. [H]) but it is still lower than that set in 
the grain mantles, again due to the formation of OH from the destruction of HDO. 

The increase in DCO"'"/HCO"'~ in the midplane between 20 and 10 AU is due to the 
increase in grain temperature which reduces the efficiency of H2, HD and D2 formation on 
grains, and therefore increases the abundance of deuterium atoms in the gas phase. The 
deuterium atoms react with HCO"*" to form DCO^. Inside of 10 AU the abundances of both 
DCO"*" and HCO"^ are very small because they are destroyed by reaction with hydrocarbons 
such as C3H4 and C2H2. 
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Fig. 4. — The fractional abundance distributions of the isotopologues of Hg in a model with 
a = 0.01. 
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Fig. 5. — The D/H ratios of the isotopologues of Hj]" in Model 1. 
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Fig. 6. — The molecular D/H ratios in the inner disk as calculated in Model 1. The D/H 
ratios are only shown in regions where the fractional abundance of the non-deuterated species 
is greater than 10"^^ (the exception is CHj]" for which D/H ratios are shown across all of 
the disk regardless of the molecular abundances) . Those species which are abundant at i? > 
10 AU, have radial scales running from 0.5 to 30 AU. Other species are only shown for the 
inner 10 AU. The atomic D/H ratio and the CH2D"'"/CHj|" ratios are important for driving 
molecular deuteration of some molecules. Relatively high molecular deuteration is retained 
even in warm regions, since the ratios are set in the cooler outer regions of the disk and 
the material is transported inwards at a faster rate than the chemical reactions that would 
destroy the deuteration. 
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Fig. 6. — cont 
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Table 5. The molecular D/H ratios in the midplane for Model 1. • ■ • indicates that the 
molecular fractional abundances are very low (< 10~^^). X:gr indicates a grain mantle 

species 



Molecule 


Input 


25 AU 


10 AU 


5 AU 


1 AU 


D/H 


2.00 


(-2) 


5.25 (-3) 


6.70 (-5) 


5.50 (-3) 


5.50 (-6) 


DCO+/HCO+ 


9.00 


(-2) 


1.89 (-2) 


2.30 (-2) 






CHsD/CH^ 


1.10 


(-1) 










H2D+/H+ 


2.00 


(-1) 










HDCO/H2CO 


7.00 


(-2) 


1.69 (-2) 


1.07 (-2) 


1.60 (-2) 




HDO/H2O 


1.00 


(-1) 








1.00 (-2) 


D2O/H2O 


1.00 


(-3) 








1.06 (-4) 


DCN/HCN 


3.20 


(-2) 








5.53 (-6) 


NH2D/NH3 


3.10 


(-2) 






4.09 (-3) 


4.50 (-3) 


HD0:gr/H20:gr 


9.00 


(-3) 


9.34 (-3) 


9.50 (-3) 


9.60 (-3) 


1.07 (-2) 


DCN:gr/HCN:gr 


9.00 


(-3) 


3.00 (-2) 


2.96 (-2) 


2.99 (-3) 




HDC0:gr/H2C0:gr 


4.80 


(-2) 


2.96 (-2) 


1.54 (-2) 


2.17 (-3) 




NH2D:gr/NH3:gr 


5.13 


(-3) 


5.34 (-3) 


5.26 (-3) 


5.26 (-3) 





- 30 - 



3.3. Ionization fraction 



Disks around young stars provide a means of dissipating angular momentum and tliey 
regulate the rate at wliicli mass accretes onto the protostar. The evolution of the star-disk 
system is controlled by angular momentum transport, but the mechanism b y which this is 
achieved is unclear. It is likely to involve turbulence (see lMorfill fc Volklll984l . and references 
therein) since molecular viscosity is too low to have much effect. The mechanism driving 
the turbulence ha s not ye t been identified b ut sev e ral processes ha v e been sugges t ed e.g . 



Balbus fc Hawlev 



Lin fc Papaloizou fllQSOl ) : IRvu fc GoodmanI fll982l): IPubruUd (19931): iBalbus et al 



Jl99lh: lLi et al.N200oh : lKlahr fc BodenheimeJ J2o3 ): lKlahij j20ol 



1996); 



( 2OO5I ): Shalybkov fc Riidigei ( 2005 ). Magneto-rotational instabil ity (MRI) is widely re 



DubruUe et al 



garded as the most promising mechanism for driving turbulence (IBalbus fc Hawleyl Il991 



Hawley fc Balbuslll99ll : iHawley et al.lll996l ). In MRI, magnetic field lines linking gas at dif- 
ferent distances from the star are stretched due to the decrease in orbital frequency with 
radius. Magnetic tension forces decelerate the inner gas, which spirals inwards and ac- 
celerates the outer gas, which spirals outwards so that the bend in the field line grows 
with time. Instability occurs if the field lines are frozen in the gas. MRI requires the 
ionization fract ion x(e) to be above a minimum value (> 10"^^ at 1 AU in the midplane 
of the MMSN Jinutsuka fc Sand (120051 )). There is some debate as to whether MRI is ap- 
plicable in disks , whe r e high densitie s and low temperatures lead to low ionization levels 
(IBlaes fc Balbuslll994l : lGammielll996l ). Many models have found that active regions (where 
the ionization level is high enough to drive MRI) can exis t in the disk along with dead zones 



where t he magnetic field is no t coup l ed to the gas e.g. ISano et al 



((20021); iMatsumura fc Pudritd (120031 ): ISemenov et all (l2004f ): llkner fc NelsonI (|2006a|Jb|ld) • 



(12000): iFromang et al. 



The size and location of the dead zone varies with the model, but most cover the region of 
terrestrial planet formation and i t is thought that dead z ones may also play a role in halting 
the inward migration of planets (IMatsumura et al.l 120071 ). 



The total ionization level in Model 1 is shown in Figure [7](a). A dead zone (defined as a 
region where x(e) < 10^^^) is present in the midplane inside of ~ 10 AU. The extent of the 
dead zone is sensitive to the mass of the disk - increasing the mass by increasing the value of 
a to 0.025 (Model 2) extends the dead zone out to 18 AU in the midplane (Fig. [7](b)). This 
is a consequence of the higher densities in Model 2, which both increase the rate at which 
ions and electrons collide and neutralize and reduce the cosmic ray ionization rate in the 
midplane. The effects of an increase in disk mass on the abundance distributions of other 
molecules are discussed in Section j 



The dominant ion varies with position in the disk. Fig [8] illustrates this for Model 1. At 
the surface, is the main ion at all radii. It is formed efficiently by X-ray ionization of H 
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atoms. Below this is a layer of which lies above a thin layer of hydrocarbon ions e.g. C4H2 
and CsHjj'. The hydrocarbon ions dominate in a region where there are still some photons 
and where the high abundance of C+ means that hydrocarbons and their ions can form 
efficiently. HCO+ is the main charge carrier below the hydrocarbon layer for R> 8 AU. It is 
destroyed at i? < 8 AU by reaction with hydrocarbon chains e.g. C3H4 and C2H2 producing 
hydrocarbon ions that are the main contributors to the charge in the central regions. HCO''" 
also contributes to the production of another abundant ion in the inner regions - CH3CO"'" 
- by reacting with CH2CO. 
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Fig. 7. — Left: (a) The ionization level in disk of Model 1. A dead zone (defined as the 
region where x{e) < 10"^^) exists in the midplane at i? < 10 AU; Right: (b) The ionization 
level in the higher mass disk of Model 2, showing an increase in size of the dead zone with 
increasing disk mass 



(a) Model 1 (b) Model 2 
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Fig. 8. — The distribution of ions in Model 1. HCO"^ is the dominant ion at i? > 8 AU, 
CH3CO"'" and hydrocarbon chain ions dominating inwards of this. At the surface the main 
ion is H"*", with a layer of C"*" below. Hydrocarbon ions exist in a layer below the C^. A 
similar distribution is seen in Model 2. 

Major charge carriers in disk 




Radius (AU) 
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3.4. Effect of disk mass 

Model 1 described above lias a = 0.01 and a total disk mass of 0.032 M© inside of 
300AU. The surface density is 38.8 g cm~^ at 5 AU. This is somewhat lower than that of the 
MMSN model and we have therefore also considered a higher mass model, again provided 
to us by Dr. P. d'Alessio. Model 2 has a higher value of a (= 0.025) resulting in a higher 
surface density (150 g cm~^ at 5 AU). 

We find that the mass of the disk does not greatly affect the molecular distributions in 
the inner disk (Fig. [9]). The higher mass does result in a narrower surface layer of molecules 
e.g. for CN, H2O etc, and the slightly higher midplane grain temperature in Model 2 means 
that molecules desorb at slightly larger radii. 

The deuteration of some molecules varies with the disk mass. Figure [TO] shows the radial 
variation of column density deuteration for both models. HDO/H2O and DCO"'"/HCO"'" 
are similar in the two models, but for DCN/HCN and HDCO/H2CO there is a drop in 
deuteration with the higher disk mass. DCN and HCN exist in a layer above the midplane 
in both models, with the Model 2 layer being at slightly higher z due to the increase in 
optical depth of this model. In both models the abundance of DCN is mainly dependent 
on the ratio of D atoms which react with HNC to form DCN, and H atoms which destroy 
DCN forming HCN. The ratio of atomic D/H is higher in Model 1 than Model 2, leading to 
a higher DCN/HCN ratio. 
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Fig. 9. — Fractional abundances in the inner disk for Model 2 (a = 0.025). The first few 
figures, from CO to CN, show those molecules which are present throughout the disk, the 
later ones show those molecules which are only present in the gas phase in the inner 10 AU. 
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Fig. 10. — The radial variation of column density deuteration (i.e. iV(XD)/A^(XH)) for 
Model 1 (dashed line) and 2 (solid line). 
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3.5. Comparison with previous inner disk models 



Markwick et al.l (120021 ) use a similar, but undeuterated, chemical network to model the 
chemistry within 10 AU of a protostar. However there are major differences between our disk 
model and theirs which result in very different temperature st ructures and hen ce different 
calculated abundances. Although our Model 1 and that used by lMarkwick et alj have similar 
grain temperatures in the midplane at 10 AU, their model has a more rapid increase in mid- 
plane temperature as the radius decreases. They also find a vertical temperature inversion 
at all radii (< 10 AU) so that the surface layers are cooler than the midplane. The relative 
surface and midplane temperatures are controlled in part by the surface density of the disk 
- the higher the surface d ensity, the thicker the disk and the warmer the midplane is relative 
to the surface layers. The lMarkwick et al.l disk has a mass accretion rate that is 10 x higher 
than ours, leading to a higher surface density and hence a colder surface temperature. 

The gra in temperature is very important since it governs the location at which molecules 
can desorb. iMarkwick et al.l find that at 10 AU almost all the molecules are accreted onto 
grains in the surface layer. Because of our higher grain temperature in this region we find 
most molecules are in the gas phase, with only water and its isotopologues remaining as 
ice. We also find that grain mantles persist at smaller radii compared to Markwick, again 
because of the difference in temperature distribution assumed. 

More recently th e chemistry in the disk within 3 AU of the star has been studied by 
Agundez et al.l (120081 ). These authors concentrated on the warm, UV irradiated, surface 
regions and found that photodissociation of CO provided a pool of carbon atoms that leads 
to the formation of simple organics such as C2H2, HCN and CH4. The formation of C2H2 
depends on the reaction of H2 with either C2 or C2li - reactions that are not generally 
included in astrochemical models since they have activation barriers of ~ 1400 K. However 
in the warm inner disk the gas temperatures are high enough for them to make a significant 
contribution to the carbon chemistry. These authors find good agreement with observations 
and their results (apart from C2H2) are in broad agreement with those presented here. 



Deuterium chemistry has previously been studied in the inner disk by lAikawa fc Herbst 



( ll999bl ) who follow the evolution of material into a radius of 30 AU. Their model differs 
from ours in several ways. They include a description of the effects of collapse on the 
input abundances. Although they include gas-grain interactions (freezeout and thermal 
desorption) they do not consider the effects of reactions on the grain surface (except for 
the formation of H2 and for recombination of ions and electro ns). This leads to a number 
of important differences between their conclusions and ours. lAikawa fc HerbstI find that 
gas phase chemistry is important for determining the deuteration of some molecules. Most 
notably this includes the DCN/HCN ratio observed in comets. They find that DCN and 
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HCN are formed by CHD + N — > DCN + C and CH2 + N — > HCN + C respectively. The 
derived ratio can be affected by both the cosmic ray ionization rate (lowering (qr results in 
a lower D/H ratio) and by whether or not HCO"^ recombines dissociatively with electrons on 
grain surfaces. In contrast we find that to a large extent the DCN/HCN ratio calculated in 
ices at the comet formation radius in our disk is driven by what happened in the molecular 
cloud, where grain surface chemistry plays an important role in determining the deuteration. 

We know from observations that the amount of water ice seen in molecular clouds can 
only be produced by grain chemistry. This alone demonstrates the importance of considering 
grain chemistry in models, but the mechanism for including it is still a matter for debate. 
Another issue that remains unresolved is the nature of non-thermal desorption processes that 
prevent the complete removal of molecules from the gas in regions where ther r nal desorption 



is inefficient. Other observ ations of star forming regions fjParise et al.l l2002l : iLoinard et al. 



2OO2I : iBacmann et al.ll2003l ) also indicate that grain chemistry is important in the formation 
of molecules such as deuterated formaldehyde and methanol and hence cannot be ignored 
when modeling deuteration effects. Equally our model falls short in not considering the 
collapse phase and the effects on the chemistry of th e incorporation of c loud m aterial into 
the protostellar disk. Neither our model, nor that of lAikawa fc HerbstI (jl999bl ) provides a 
complete picture of the chemical history of disks. 
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4. Comparison to observations 
4.1. Cometary ices 



The comets we see today are the remains of icy planetesimals that formed early in 
the history of the solar system. Since their formation they have been kept mainly at cold 
temperatures which restrict further chemical evolution. Consequently their chemical make 
up reflects to first order, the composition of the protostellar disk at their time of formation. 
Simple species such as H2O, CO, CO2, HCN, CH3OH, H2CO have been observed in the 
comae of many comets, but D/H ratios have been measured in only four comets - Halley 
(HDO/H2O), Hyakutake (HDO/H2O), Hale-Bopp (DCN/HCN, HDO/H2O) and C/2002 T7 
(LINEAR). These are all long-period comets, and so formed inside of the Trans- Neptunian 
region, after which their orbits were perturbed by interaction with the giant planets and they 
were ejected outwards towards the Oort cloud. The chemical composition of the comets 
should reflect the temperature, density and ionization state of the region in which they 
formed. Measurements of the D/H ratios of cometary ices could also be used to constrain 
the degree of processing of the interstellar ices incorporated into the protosolar disk, and 
place limits on the amount of radial mixing experienced by material in the disk. 

The observed molec ular deuteration is similar for all four comets. In Hale-Bopp, 
DCN/ HCN = 2.3 xlQ-^ jMeier et al.lllQQsh . (D/B. )ur) is measured as 3.3 xlO"^ in Hale- 



Bopp (IMeier et al. 



1998 



xlO"'^ f Eberhardt et al. 



) 2.9 xlQ-'^ in Hyakutake flBockelee-Morvan et all 119981 ) and ~ 3.1 



I995I: iBalsiger et al.lll995h in Co met/P Halley In C/2002 T7 (LIN- 



EAR) (D/H)jj^Q is 2.5 xlO ^ ( Hutsemekers et alll2008 ). (These measurements correspond 
to x(HDO) /x(H20) ratios of 6.6 x 10-^ 5.8 x 10-^ 6.2 x 10"^ and 5.0 x 10"^ for Hale-Bopp, 
Hyakutake, Halley and C/2002 T7 (LINEAR) respectively) 

Figure [TT] shows the variation in the x(DCN)/x(HCN) and x(HD0)/x(H20) ratios in the 
ice mantles over the radius range from 0.5 to 10 AU. Our model shows very little variation in 
the ratios in the midplane in this region. The x(HD0)/x(H20) ratio is set in the molecular 
cloud phase and very little processing of the water ice occurs in the disk. Hence in the 
midplane, the x(HD0)/x(H20) ratio in ices has a constant value of 0.022 from 30 to 2 AU. 
DCN/HCN shows some slight variation but the value is still roughly 0.018 from 5-10 AU. 
For both molecules the model deuteration ratios are much higher than observed. The model 
DCN/HCN is approximately 8 times higher than observed, and the HDO/H2O ratio is 33 - 
44 times higher than observed. Hence, although comets have to contain ices accumulated at 
low temperatures they must also include some ices that were formed in a higher temperature 
phase. 



There are several possible explanations for the discrepancy between our model and the 
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observations. Since the disk chemistry does not change the HDO/H2O ratio from the value 
set in the molecular cloud it may be that the problem lies in our molecular cloud model, and 
that HDO ice is formed too efficiently. In our cloud model the ratio of HDO/H^ O in the ice is 
0.9 %, compared to an upper limit of ~ 0.5 - 2 % towards low mass protostars fjParise et al. 



2OO3I ). This upper limit is much less than the deuteration ratios seen in methanol and 
formaldehyde which are also expected to form on grain surfaces. If all three molecules form 
at the same time then their D/H ratios should be similar and should reflect the gas phase 
atomic D/H ratio at the time of formation. Observations of methanol and formaldehyde in 
low mass star formation regions find much higher D/H ratios that the upper limit derived 
for HDO/H2O ice, suggesting that there are still unresolved issues in our understanding of 
the formation of these molecules. 

Another possibility is that the ices are removed from the grains during infall through 
the accretion shock at the surface of the disk. Models suggest that most of the water ice on 
infalling grains can be remo ved within 30 AU of a solar ty pe star (e.g. 



Lunine et al. 1991 



Neufeld &: Hollenbachlll994l ). Recent Spitzer observations (IWatson et al.l 120071 ) have found 
water vapor emission in a Class object (NGC 1333-IRAS 4B), which has been attributed to 
emission from a surface accretion shock. Water would therefore enter the disk as vapor rather 
than as ice. This could affect the D/H ratio of the ices that reform in the disk, since water 
formed in the warm post-shock region would have a lower D/H ratio than water formed in 
the cold molecular cloud and this ratio would be transferred to the ices when the molecules 
were redeposited. However it is unclear how widespread or long-lived this phenomenon is 
since water emission was detected in only one object out of a sample of 30. 

Although vertically av eraged hydrod y namical disk mod e ls fin d that material is ad- 
vected radially inwards e.g. iPringld (jl98ll ): iLin fc Papaloizoul (119851 ). with outward fiow oc- 
curring only in the cold outer regions of the disk, more complex models have shown that 



Fig. 11. — The deuteration of HCN and H2O ices in Model 2. 
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large scale meridional flow patterns develop in d i sks e.g. lUrpinI (1198 



Rozvczka et all fll994f ): lRegev fc GitelmanI (120021 ) : lTscliarnuter fc Gaill (120071 ). In these mod 



\: iKlev fc Linl (Il992h : 



els material flows inwards at the surface and outwards in the midplane, leading to a mixing 
of warm and cold material. Observational evidence of the importance of mixing in disks 
comes from comets. Silicates are thought to have undergone signiflcant high temperature 
processing, which converted amo rphous silicate s into the crystalline form, before they were 
incorporated into comet nuclei (IHannerl Il999l ). One way to do this would be by radial 
mixing where the silicates were heated to high temperatures in the inner disk before being 
mixed outwards to the comet formation region. Additi onal evidence of mixing in the history 



of co mets has been seen in the STARDUST samples (IBrownlee et al.ll2006l : IZolensky et al. 



20061 ). If grains can be mixed then gases could be too, and this would have the effect of 
changing the deuteration ratios by bringing together molecules formed in the cold outer disk 
or in the molecular cloud with those formed in warmer regions of the disk. Overall, mixing 
could be expected to reduce the D/H ratios seen in cometary ices. 

Diff usion has been includ ed i n models of the deut e rium chemistry in th e midplane of 
disks by iDrouart et al.l (119991 ) and iMousis et al.l (120001 ). iDrouart et al.l (119991 ) assume that 
the water is initially gaseous with (D/H)jj q = (D/H)jj^ and t hey follow the evolu tion of 
(D/H)]^ Q over time in a model that includes turbulent diffusion. iMousis et al.l (l2000l ) take a 
similar approach and modeled the DCN/HCN ratio. These authors concluded that cometary 
ices condense in the nebula, rather than being retained during the disk and star formation , 
and are composed of material that is partially reprocessed in the disk. iHorner et al.l (120071 ) 
also found a relationship between the D/H ratios in ices and the distance from the Sun at 
which they condense. 

Our model is therefore probably too simplistic in its assumptions and further work is 
required to study how the accretion shock alters the D/H ratios and the survivability of th e 
ices in the star forma tion p rocess, and to consid er the effects of mixing. Illgner et al.l (120041 ): 
Semenov et al.l (120061 ) and IWillacy et al.l (120061 ) have all shown that mixing in the vertical 
direction can greatly affect the abundances in the disk, but isotop e s have not ye t been 
stu died with these mod els. Also, as shown by IDrouart et al.l (119991 ): IMousis et al.l (l2000l ) 
and lHorner et al.l (120071 ). radial mixing alters D/H ratios in ices, bringing cold material into 
the inner disk where it can desorb and returning material to the colder regions where it can 
reaccrete. Both radial and vertical diffusion alter the molecular D/H ratios by combining 
material formed under different temperature regimes. 
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4.2. Gas phase observations of extrasolar disks 



Our basic model (excluding isotopes) is the same as used in IWoods fc Willacyl (120091 ) 
who considered the chemistry of carbon isotopes. That paper includes a comparison of the 
result to recent infrared observations of the inner regions of protostellar dis ks. Here we merely 
summarize our findings for completeness, but the reader is referred to IWoods fc Willacy 
(120091 ) for more details. 



Several molecules e.g . CO, C2H2, HCN, OH, H^O and COo have been obse rved in AA 
Tau JCarr fc Naiitalboosi), PR Tau an d AS 205 A JSalvk et allboosh . GV Tau Joibb et al. 



20041 . 12OO7I ) and IRS 46 (ILahuis et al.ll2006l ). The derived rotational temperatures suggest 
that the molecular emission comes from a region above the midplane. In our models, these 
temperatures are reached in the region where thermal desorption is efficient, and the gas 
and grain temperatures are decoupled. We find that in Model 1, CO, H2O and OH have 
peak abundances in regions of the disk with similar temperatures to the observed rotation 
temperatures for these molecules. The surface peak of CO2 also coincides with the observed 
rotation temperature, although the peak abundance of this molecule is closer to the midplane. 
For C2H2 the observations lie above the modeled distribution, and for HCN they are at the 
upper edge of the model. 

We can also compare our results with the column densities observed (Table E]). The 
calculated column densities at 1 AU for several molecules are given in Table H] and are much 
higher than observed, but these values are calculated from the disk midplane to its surface 
and therefore include material that is much colder than observed. If we instead calculate the 
column densities only in the upper regions of the disk where the gas temperature is warmer 
than 150 K, the agreement is better. The calculated column densities are N(HCN) = 4.4 
xlO^^ cm-2, N(CO) = 1.6 xlO^^ cm'^, and N(C02) = 1.3 xlO^^cm'^ for Model 1. N(C2H2) 
is far too low (9 xlO® cm~^). This low abundance may be due to the exclusion of reactions 
involvin g H^ with C2 or C^H that were found to be important in the formation of C2H2 in 
disks by lAgundez et al.l (120081 ) but are not included in our reaction scheme. The calculated 
N(H20) is also too high compared to observations, but N(OH) gives good agreement with 
the abundances measured in DR Tau and AS 205A ISalyk et al.l (120081 ). 



A further test of the models is to compare the C H4/CO ratio. iGibb et al.l (120071 ) found 
an upper limit of 0.0035, similar to the upper limit iGibb et al.l (120041 ) found for this ratio 
in HL Tau (< 0.005). The calculated ratio at 1 AU of 0.006 is consistent with both these 
observational measurements. 



Oi et al.l (I2OO8I ) observed deuterated molecules in TW Hya and found a gradient in the 
DCO+/HCO+ ratio. The ratio increases from 0.01 at 30 AU to 0.1 at ~ 70 AU. There is a 
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steep fall-off in N(DCO"'") at i? > 90 AU. Our current r aodel only ext ends out to 35 AU but 
when we consider it together with our outer disk model (IWillacyll2007l ) we can see if the same 
trend shows up in the chemical modeling data. Table [7] shows the model column densities 
for DCO"*" and HCO+, as well as their ratio for radii between 30 and 150 AU. Although we 
do not find the same steep fall off in N(DCO''~) at i? > 90 AU, we do see a similar increase 
in DCO"'"/HCO"'" from 30 - 100 AU followed by a fall at larger radii. However our D/H ratio 
at 100 AU is ~ 10 xlarger than inferred from the observational data. 
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Table 6. The observed column densities in the inner regions of protostellar disks. The 
values for Model 1 are calculated at lAU for the surface layers of the disk where T 



gas 



> 



150 K. 



Molecule 




Column Density (10^^ cm" 








GV Tau^ 


IRS462 


AA Tau^ 


AS 205A4 


DR Tau^ 


Model 1 


H2O 






65.0 


60.0 


80.0 


1280.0 


OH 






8.1 


20.0 


20.0 


23.0 


HCN 


3.7 


5.0 


6.5 






4.4 


C2H2 


7.3 


3.0 


0.8 






9.0 (-8) 


CO2 




10.0 


0.2 






13.0 


CO 


590.0 


200.0 


49.0 


60.0 


70.0 


1600.0 



References. 



Gibb et al. ( 



Carr fc Naiital ((20081) (4) ISalvk et all fl2008f ) 



200'i 



3), (2) iLahuis et all J2006h . (3) 



Table 7. The calculated values of N(DCO+), N(' HCO+) and DC O+/HCO" 

i? > 30 AU is taken from Iwillacvl J2007h . 



The data for 





30 AU 


50 AU 


100 AU 


150 AU 


N(DCO+) 


1.9 (12) 


1.5 (12) 


6.6 (12) 


2.3 (12) 


N(HCO+) 


3.3 (13) 


1.2 (13) 


3.9 (12) 


2.8 (12) 


DCO+/HCO+ 


0.05 


0.125 


1.7 


0.8 
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Conclusions 



We have presented the results of calculations of the chemical abundances in the inner 
regions of protostellar disks including deuterium chemistry. We find that for many molecules 
the deuteration is set by cold temperature chemistry in the parent molecular cloud, and that 
the chemistry in the disk itself has little effect. In particular the calculated D/H ratios in ices 
in the region of cometary formation are found to be high relative to the (limited) cometary 
observations. This suggests that our model is not complete, and additional processing of 
ices is required after the molecular cloud stage, beyond the chemist ry included in our model s 
here. Our carbon isotope modeling leads to a similar conclusion (jWoods fc Willacyll2009l ). 
The processing could be achieved in the accretion shock, and/or by mixing in the disk, both 
of which could lead to lower D/H ratios. Additional work is required to fully understand 
the process by which interstellar grains and ices are incorporated into planets and other 
planetary bodies. 



When a comparison is made with the observations of external disks (jCarr fc Najita 



2008; 


Gibb et al. 


2007. 


2004) 



in our disks and that derived from the observations. The molecular emission comes from a 
region above the midplane, where thermal desorption is efficient. 

We find major differences in the results of our calculations and those of other models, 
mainly as a consequence of either the physical disk model assumed, or because of the chemical 
processes include d . We find that the disk chemistry has less effect on D/H ratios than 
Aikawa fc HerbstI (jl999bl ) because we allow chemistry to occur on the surfaces of dust grains 
and hence gas phase chemistry is less important in determini ng deuteration. The tem perature 
structure in our models is very different to that assumed by lMarkwick et al.l fl2002l ). so that, 
despite using similar chemical networks, we get very different abundance distributions in 
the inner 10 AU of the disk. This highlights the effect that disk structure can have on the 
molecular distributions in disks. As observations improve in resolution, for example with the 
advent of ALMA this may enable the differences between models to be used to distinguish 
among them observationally. 
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